SMALL POLYGONS AND TORIC CODES 



GAVIN BROWN AND ALEXANDER M. KASPRZYK 

Abstract. We describe two different approaches to making systematic classifications of plane 
lattice polygons, and recover the toric codes they generate, over small fields, where these match 
or exceed the best known minimum distance. This includes a [36, 19, 12]-code over F7 whose 
minimum distance 12 exceeds that of all previously known codes. 



1. Introduction 

We are interested in planar lattice polygons P - not least for their own sake - and, for a finite 
field Fq, the toric codes Cp{¥q) that they determine. We recall the construction of such codes, 
following Hansen |Han001 IHanl2j . in Section [2] below. Briefly, P and q together determine a 
k X n integer matrix - the toric code - for integers k, the dimension, and n, the block length. 
The important invariant of a code C is its minimum distance d, the shortest Hamming distance 
between any two distinct points of the image of the matrix - that is, the distance between points 
of the lattice generated by the rows of the matrix. It is this measure which limits the ability of 
the code to detect and correct errors. 

The collection of all toric codes includes some champions among linear codes, in the sense 
that they have minimum distance greater than that of any other known code with equal block 
length and dimension. This is exemplified by the paper of Joyner |Joy04| and the recent paper 
of Little [Litllj . the latter describing a generalised champion [49, 12, 28]-code over Fg (we follow 
the standard notation and denote the invariants of a code by a triple [n, k,d]). We find a new 
champion, a [36, 19, 12]-code over F7, however the novel part of our paper is the systematic 
approach to classification of polygons that we use to discover this example. This approach 
exploits database methods in collaboration with computer-aided polygon computations in a 
way that seems not to have been pursued in this area. As a spin-off, we classify a certain set of 
'small' lattice polygons, as we now explain. 

One of the key points when considering toric codes over Fg is to restrict attention to polygons 
P that lie in a square of side-length m < q — 2. In this paper we concentrate on fields Fg with 
q < 9 a prime power, and so we first enumerate all polygons that lie in a square of side-length 
m < 7. We regard two polygons P and P' as equivalent if there exists a point it G and change 
of basis if of the underlying lattice such that ip{P — u) = P' . Table [1] summarises our results by 
listing the number of integral polygons (up to equivalence) that lie in the square [0,m] x [0,m], 
but not in a smaller square. It also lists the maximum number of vertices that such a polygon 
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can have, and the number of polygons that achieve this upper bound. The complete hst of 
polygons is available on the webpage |BKj . 

The systematic approach described above would certainly exhaust all toric codes over small 
fields (even the generalised codes of |Rua091 ILitllj ). but it obviously becomes less practical as 
m increases. It is remarkable to us that the number of polygons grows as slowly as it does, 
although this is presumably an artefact of two dimensions. We seek champion codes, that 
is, codes exceeding the current largest known minimum distances. Computing the minimum 
distance is very often too hard (it is certainly computable, but in many cases we have reliable 
estimates of computation time of the order of a year), but determining whether or not a given 
code is a champion is much easier: in most cases, simple linear combinations of the rows of 
the corresponding matrix provide sufficiently short vectors to rule it out. We carry out this 
calculation for each one of these polygon codes for every applicable field of size at most nine. In 
this early part of the polygonal classification, there is exactly one champion - the code mentioned 
above - as well as many examples that match the current best known minimum distance. All 
but 257 of the polygonal codes have a short vector that is a combination of at most three rows of 
the generator matrix; only six of these, other than the champion, require more than four rows. 

Table 1. The number of polygons (up to equivalence) contained in the square 
[0, m] X [0,m], but not contained in a smaller square. 
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^Polygons 
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7^ Max polygons 
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1 


1 


1 


15 
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The champion polygon and many of the equal champions are LDP-polygons (lattice polygons 
associated to log del Pezzo surfaces, explained in detail in Section H]). Motivated by this, we 
go on to consider codes associated to toric log del Pezzo surfaces, using the classification of all 
15,346 LDP-polygons of index i < 17 [KKNlOj . The difficulty here is that these polygons are not 
presented as lying in the smallest possible square, so we describe an algorithm that determines 
the smallest square that contains a polygon equivalent to a given polygon. This algorithm - an 
adaptation of the minimum width algorithm - is the main point, since our search through the 
resulting data for toric log del Pezzo codes over the smallest possible field reveals no further 
champions. 

2. Toric codes 

Let M = be a two-dimensional lattice, and let P C M (S)z Q be a convex lattice polygon 
with vertices in M. Let ¥q be a finite field such that, up to translation, P is contained in the 
square conv{(0, 0), {0,q — 2), {q — 2,0), {q — 2,q — 2)}. Hansen demonstrated in jHanOOl rHanl2j 
how one can associate a linear code to P. Let e € Fg be a primitive element of the field, and for 
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each lattice point u := (ui, U2) E P PI M define: 

e{u) : F* X ^ Fg 

The set of vectors 

{ie{u){e\e^))o<i,,<g^2 \ u e P D M} 

then generates a Unear code of block length n = {q — 1)^ and dimension k = \P n M\, denoted 

by Cp{¥g). 

We sketch the connection with toric geometry. For the details see |Han001 IHan021 Joy04| . 
Associated with any polygon P is a non-singular fan A in := Hom (M, Z) given by the 
following construction. Let hp : Nq — )• Q be the support function defined by 

hp{v) := inf{(n, | u S P}. 

This is a piecewise linear function, and partitions Nq into a finite collection of strictly convex 
polyhedral cones; in other words, hp defines a fan A' in Nq. We can refine A' to a non-singular 
fan A by inserting rays via a well-known process (see, for example, |Oda78] ). Notice that, for 
any cone o" G A, there exists a lattice point /q- G M such that 

hp{v) = (/ct, v) , for all v G a. 

Let Xp be the complete smooth toric surface associated with A, where the algebraic torus 
T/v is defined by T/v := Hom ^M, F*^ . The polygon P is associated with the Cartier divisor 

Dh-=- ^ (/p,/)) • vert(/>) , 

Rays p of A 

where vert(/9) = is the closure of the T/v orbit of p in Xp (again, see |Oda78j for the details) 
as follows. The space of global sections H^{Xp,OxpiDh)) has dimension |PnM| with basis 
{e{u) \ u G P n M}. Thus the code Cp{¥g) can be obtained via: 

The generators of the code are given by the image of the basis; i.e. 

{(e(n)(t))jgT^ \ uePnM}. 

The toric codes Cp(Fg) obtained via translation of P, or more generally by any affine linear 
isomorphism of M = Z^, are monomially equivalent |LS071 Theorem 4], hence it is sufficient 
to consider P up to equivalence. Following jRuaOQl ILitllj one can restrict the evaluation map 
(j2.ip to a subspace of the Riemann-Roch space to construct generalised toric codes. 



4 



G. BROWN AND A. M. KASPRZYK 



3. Classifying small polygons 

Definition 3.1. Let P C A/q he a lattice polytope, and let v S vert(P) be a vertex of P. If 
the polytope P^ := conv((i-' fl M) \ {v}) satisfies dimP^ = dimP, then Py is said to have been 
obtained from P by shaving. Given a polytope Q we say that Q can be obtained from P via 
successive shaving if there exists a sequence of shavings Q = P^'^^ C • • • c = P, where 
p(i-i) ^ pW g^^^ ^. ^ vert(p(*)), <i<n. 

Lemma 3.2. Any lattice polygon Q in the box Bm := [0,m] x [0, ?7i] can be obtained from Bm 
via successive shaving. 

Proof. Let Q C Bm- If Q 7^ Bm, then we can certainly shave a vertex v of Bm to obtain a 
polygon Py that contains Q. Continuing inductively, suppose P C Bm is a polygon with Q C P. 
li P ^ Q then there exists some vertex v € vert(P) such that v ^ Q and Q C. P' := Py But 

Vol(Q) < Vol(P') < Vol(P) . 

Since these normalised volumes are integers, this process must terminate after finitely many 
steps. □ 

Remark 3.3. The normalised volume Vol(P) of a lattice polygon P is twice the Euclidean volume 
vol(P); when P is a lattice polygon, vol(P) € (1/2)Z. 

It is clear how this lemma presents an algorithm for classifying small polygons. Choose a box 
Bm and initialise a sequence containing that box as a polygon: 

Ps: = [Polytope ([[0,0] , [0,m] , [m,0] , [m,m]])] ; 

The algorithm now extends Ps recursively by shaving the vertices from each polygon in the 
sequence, checking for equivalence before appending each new polygon. 

idx:=l ; 

while idx le #Ps do 

for V in Vertices (Ps [idx] ) do 

P : =Polytope (Exclude (Points (Ps [idx] ) , v) ) ; 

if Dimension (P) eq 2 and not &or [IsEquivalent(P,Q) : Q in Ps] then 

Append (~Ps,P) ; 
end if ; 
end for; 
idx +:= 1; 
end while; 

The appeal to IsEquivalent is simply checking that the new polygon P is not equivalent to a 
polygon already contained in the sequence. This is the basic algorithm (after checking simple 
invariants, attempt to build an isomorphism working up from low- valency vertices), which works 
well for small values of m; larger boxes require some straightforward optimisations that we do 
not describe. 
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4. Log del Pezzo codes 



A polygon P C N (8>z Q is said to be an LDP-polygon if it contains the origin strictly in its 
interior, and if each vertex of P is a primitive lattice point in N := Horn (M, Z). In this case the 
corresponding toric variety is a log del Pezzo surface: it has only cyclic quotient singularities 
and —Kx is ample. For example, the polygon conv{(l,0), (0, 1), (—1, —1)} is an LDP-polygon 
corresponding to the projective plane P^. For such P, the dual polygon P* is defined by 



Typically P* is not a lattice polygon, but there exists a smallest integer £ > 1 for which the 
vertices of iP* are lattice points. This integer i is called the index of P. It equals the Gorenstein 
index of the corresponding log del Pezzo surface. 

There are infinitely many LDP-polygons, but finitely many for any fixed index i. Figure [T] 
shows an LDP-polygon with vertices (—1,2), (1,2), (1,0), (—1,-1), (—2,-1). It has index 
i = 10; the dual is a lattice polygon after dilating by a factor of 10. 



Figure 1 . An LDP-polygon of index i = 10, with its dual rational polygon. 

The classification of LDP-polygons of indices 1 < i < 17 is given in [KKNlOj and available 
online at [BK] : the example in Figure [1] is No. 523 in the classification. To use these for 
toric codes, we allow for translation, and must find equivalent polygons that lie in a box Bm, 
rather than containing the origin as given. The representative polygons in the database are not 
necessarily the smallest, so we describe an algorithm to determine, for given polygon P, the 
smallest m for which there is a polygon P' C Bm equivalent to P. (We treat LDP-polygons 
simply as a source of interesting examples to generate codes, and therefore regard P as lying in 
Mq, rather than in terms of algebraic geometry, where they should be viewed as lying in Nq. 
Our work here is entirely in terms of Mq, so no confusion arises from the name.) 

Lemma 4.1. There is a constructive algorithm that, for any given polygon P C Mq, determines 
the smallest m for which there is a polygon P' C B^ equivalent to P , without enumerating any 
other polygons. 

Proof. We use a modification of the standard minimal width algorithm. Begin by considering 
Wi := wdg* {P) for each dual basis vector e* N := Hom (M, Z). li Wi < m for all i, then, after 
possible translation, P is contained in Bm. 

Suppose without loss of generality that wi > m. Translate P so that the (rational) point 



P* := {v G Mq I {v, u) > -1 for ah u e P} . 





6 



G. BROWN AND A. M. KASPRZYK 



is at the origin. Let O C P be the largest disc, centred at the origin, that is contained in P. 
Then d := wdg* (O) satisfies < d < wi. In the dual space Nq, consider the disc D centred at 
the origin of radius wi/d. We claim that wd^ (P) > wi for any primitive lattice point u £ N\D. 
Indeed 

wd„ (P) > wdu (O) > wdu' (O) = — ■ d = wi, 

where u' is the unique rational point lying on the boundary of D in the direction of u. Thus P 
is equivalent to some P' C if and only if there exists a dual basis Bm contained in D such 
that wdy (P) < m for each v € Bm, and this is a finite search. □ 

Applying this to the LDP-polygons, we find they distribute over m as in Tabled we have cut 
the table short, but it could be continued to list all 15,346 polygons, with maximum m = 68. 

Table 2. The number of LDP-polygons (up to equivalence) of index i < 17 
contained in the box [0, m] x [0, m], but not contained in a smaller box, for small 
values of m. 



m 


2 


3 


4 


5 


6 


7 


8 9 


10 


# 


11 


62 


364 


591 


1125 


777 


1277 904 


1187 



5. The computations and champion polygons 

We use the computer algebra system Magma |BCP97] . and in particular its convex poly- 
tope [BBKllj and linear codes packages. The polygons are harvested from the Graded Ring 
Database |BKj . which can be queried from within Magma via an XML interface. Since the 
number of calculations are rather large, they were run in parallel on a cluster of 144 processor 
cores at Imperial College, London. 

The LDP-polygon number 75 in the Graded Ring Database gives a code over F7 that has 
length 36, dimension 19 and minimum distance 12 (see Figure [21 this corresponds to a unique 
LDP-polygon, since the only point that could serve as an origin is (2,2)). The previous best 
known minimum distance was 11, according to Grassl's online database of linear codes |Gra] 
(also accessible from within Magma); the theoretical maximum is 15. 




Figure 2. The translated LDP-polygon with database id 75. 
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A curiosity about polygons falls out of the classification. Among those having the maximal 
number of vertices for their box size, there is a unique one of minimal volume that is homogeneous 
(by which we mean that all vertices are locally isomorphic to one another). We sketch the 
polygons with maximal number of vertices for m < 7 in Figure O notice that for m = 7 there 
are two polygons, but only the first is homogeneous. 




Figure 3. The polygons with maximal number of vertices for each box size 
which have the smallest volume. 

We reinterpret this geometrically as the following statement, which holds for box size up to 
seven. 

Conjecture 5.1. Fix m € Z>o. Consider all toric surfaces S with —Ks ample (hut no con- 
ditions on the singularities) for which the anticanonical polygon of monomials in H^(S,—Ks) 
lies in an m X m box. Among such S, consider those with maximum Picard rank. Then among 
these there is a unique surface S which is both smooth and has minimum anticanonical degree. 

6. Higher prime powers 

We cannot calculate minimum distances, or indeed champions, for higher prime powers: for 
many polygons with the potential for a large minimum distance, the computation time required 
to prove this is in the order of millions of years. For specific polygons there exist results on the 
minimum distance for general prime powers q |LS07l IYZ09j . and more can be said when the 
polygon is Minkowski decomposable |LS06j . however we do not pursue these techniques here. 

Instead, we begin to estimate minimum distances using linear combinations of four rows of 
the generator matrix. Recall from the introduction that amongst the 1,249,439 toric codes 
generated by polygons in boxes of size up to m = 7, all but six are either shown not to be an 
equal champion, or to have a prohibitively short vector witnessed by some linear combination 
of rows involving at most only four rows. 
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We work with the next two possible values q = 11 and 13, with box sizes 8-9 and 10-11 
respectively. In each case, we find all LDP-polygons from the database (of index < 17) whose 
minimum box has size m for which q is minimal satisfying q — 2 > m. In other words, each 
known LDP-polygon will be considered in a minimum box and with respect to its minimum q. 

We partition these LDP-polygons by their number of points, k. For each /c, we work through 
the polygons and find the length of the shortest vector that is a linear combination of up to four 
rows of the generator matrix, obtaining an upper bound d\y on the best minimum distance that 
could be achieved. These bounds are listed in Tables [SHU 

Remark 6.1. In theory the bound d-^ depends on the choice of embedding of P in the box Bm] 
some embeddings may give sharper bounds than others. Lemma 14.11 can generate all possible 
embeddings in B^, and experimentation suggests that for practical purposes fib is insensitive to 
the choices made. 

Table 3. The number of LDP-polygonal codes with (minimal) q = 11. The first 
row is the dimension k of the code; the second, the length of a known vector, and 
so an upper bound for the largest minimum distance; the third, the number of 
polygons (up to equivalence) having such a vector achieving this bound. 
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